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Abstract 

Hydro dynamic equations for an inelastic Maxwell model are derived from 

the inelastic Boltzmann equation based on a systematic Chapman-Enskog 

perturbative scheme. Transport coefficients appear in Navier-Stokes order 

have been determined as a function of the restitution coefficient e, which 

cannot be defined for small e as a result of the high energy tail of the velocity 

distribution function obeying a power law. The dispersion relations for the 

linearized equation around a homogeneous cooling state have been obtained. 
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I. INTRODUCTION 



The history of statistical mechanics begins with the gas kinetic theory of elastically 
interacting particles. It is well known that the Maxwell-Boltzmann distribution for the 
velocity distribution function (VDF) plays a key role in statistical mechanics [1], and the 
relaxation of nonuniform dilute gases to the equilibrium state is described by the Boltzmann 
equation [2] . Hydrodynamic equations have also been derived from the Boltzmann equation 
based on the Chapman- Enskog method [3]. On the other hand, when there are inelastic 
interactions among particles, the behavior of collections of particles is completely different 
from that of elastic particles: There are no equilibrium states and any spatially homogeneous 
states are no longer achieved. Such a collection of the particles having inelastic interactions 
is called the granular gas whose physical realization can be observed in rings of planets, 
small planets, suspended particles in fluidized beds, and rapid granular flows etc. [4]. 

A proper set of hydrodynamic equations for granular fluids depend on situations. In 
some cases, granular particles are condensate and the fluid obeys Bagnolds' scaling in which 
shear stress is proportional to the square of strains. This flow is called a 'frictional flow' and 
achieved in many situations under the gravity [5]. The derivation of Bagnolds' scaling is 
not difficult if we adopt a phenomenology [6] but is complicated if we wish to derive it from 
a microscopic motion of particles. On the other hand, there are some ffows dominated by 
binary collisions among particles. We call this flow a 'coUisional flow' or a 'rapid granular 
flow' and its constitutive equation is similar to that of the Navier-Stokes equation, i.e., the 
shear stress is proportional to the strains [7]. For rough particles near the boundaries, we 
may need to consider the effects of spins of particles [8] , but for smooth particles we do not 
need such the complications. 

For coUisional granular ffows, we may have a standard procedure to derive hydrodynamic 
equations starting from the inelastic Boltzmann equation. The most of papers assume in- 
elastic hard-core collisions among particles. In fact, there are some important investigations 
along this line. Jenkins and Savage [9] assume that the velocity distribution function obeys 
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the local Maxwellian, and derive a set of hydrodynamic equations. Later, Lun et al. [10], 
and Jenkins and Richman [1 1] remove the ansatz of the local Maxwellian and derive hydro- 
dynamic equations. Recently, Sela and Goldhirsch [12] indicate the insufficiency of the Grad 
expansion by Jenkins and Richman [11], and they have developed a systematic expansion 
of spatial inhomogeneity and small inelasticity. In practical sense, the method developed by 
Brey et al. [13] and Garzo and Dufty [14] is the most useful, in which they have extended the 
Chapman- Enskog method to gas particles with inelastic hard-core coUisions. Their method 
seems to be reasonable and can be used for any inelasticities. However, the method adopted 
there contain some unclear points; (i) All of methods contain two types of perturbations. 
One is for spatially inhomogeneity and another is for inelasticities. Although the result 
by Brey et al. [13] is comparable with numerical results with any inelasticities, the logical 
support of their method for inelastic particles far from the elastic limit is unclear. In fact, 
they assume that the VDF in a homogeneous state obeys a function with the lowest order 
Sonine expansion, which should be only valid for nearly elastic cases, (ii) Their analysis 
does not contain any information of tails in VDF which is believed to obey an exponential 
function and cannot be described by the expansion by the Sonine polynomials. The effect 
of tails may be small but it is unclear how it affects hydrodynamics. 

The difficulties in the analysis of inelastic hard-core particles come from the form of 
the coUisional integral which is proportional to the relative speed of colliding two particles. 
It is known, however, that the coUisional integral of the Maxwell molecules which has the 
potential r~*^ with the relative distance r of colliding two particles is independent of the 
relative speed [1,3,15]. Thus, the Maxwell model has been used for the analysis of the 
Boltzmann equation as the simplest model. Quite recently, this model is applied to kinetics 
of inelastic particles [16]. It is remarkable that a tail of the scaled velocity distribution 
function in a homogeneous cooling state of the inelastic Maxwell model can be obtained 
analytically [17-20]. The result of their analysis in which the tail of VDF obeys a power law 
is also interesting, because this result is consistent with the result of simulation in fluidized 
beds which has a long-ranged hydrodynamic interactions [21], and leads to some singularities 
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in higher order moments. 

The objective of this paper is to derive a set of hydrodynamic equations from the inelas- 
tic Boltzmann equation systematically. For this purpose, we adopt the inelastic Maxwell 
model, and the Chapman-Enskog method developed by Brey et al. [13]. The organization 
of this paper is as follow. In the next section , we explain the derivation of hydrodynamics 
from the inelastic Boltzmann equation for the inelastic Maxwell model. The result contains 
the transport coefficients such as the viscosity, the thermal conductivity, and the 'diffusion' 
coefficient as a proportional constant to the density gradient in the heat flux. The detailed 
calculation for the homogeneous solution and the framework of the Chapman-Enskog solu- 
tion, and explicit calculation of the transport coefficients are presented in Appendices A, 
B, and C, respectively. In Section III, we discuss the linear stabihty of the homogeneous 
cooling state. This requires consideration of the explicit form of VDF at the first order and 
some transport coefficients in Burnett order. Appendices D and E provide the calculations 
of these quantities. In Section IV, we discuss our result. In particular, we demonstrate 
that the scaling solution in a cooling state is equivalent to a steady state with the Gaussian 
thermostat. So the generalization of the analysis developed in this paper and others is rele- 
vant to discuss a steady systems of particles suspended by fluid flows. In Section V, we will 
conclude our results. 

II. BOLTZMANN EQUATION AND CHAPMAN-ENSKOG SOLUTION 

A. Framework of the Chapman-Enskog method 

In this section, we derive the Chapman-Enskog solution of inelastic Boltzmann equation 
with inelastic Maxwell particles. Here we adopt the Chapman-Enskog method developed by 
Brey et al. [13] which is the most useful in the practical sense. We also restrict our interest 
to the case of three dimensional dilute gases. Thus, the argument is parallel to that by Brey 
et al. [13]. 
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Let us assume low density gases of smooth identical particles with the mass m, the veloc- 
ity V and the diameter a. Hydrodynamic variables to characterize the macroscopic behavior 
of the gas are the number density n(r, t), the velocity field u(r, t) and the temperature field 
r(r, t) defined by 

n{r,t)^ ldvfir,v,t), (1) 
n(r,t)u(r,t) = J rfvv/(r,v,t), (2) 
^n{r, t)T{r, t) ^ J dv^mV'f{r, v, t), (3) 

where V = v — u and /(r, v, t) is the position and velocity distribution function. 

The distribution function /(r, v, t) in our system obeys the inelastic Boltzmann equation 

(9, + vV)/ = J[/,/]. (4) 

The coUisional integral in the inelastic Maxwell model becomes 

J[f,h] = a\J'^ [ rfvi / d&{e-'b-' - l)f{r,v,t)h{r,v^,t) (5) 

where Tq is the characteristic temperature for the potential among particles, g = v — Vi, 
a is the unit vector along the line connecting centers of mass of contacting particles. The 
operator is the inverse of the coUisional operator h which are defined as 

^'g = g- (l + e)(g- <t)«t (6) 

r^g = g-ii^(g-<T)(T, (7) 

where e is the coefficient of restitution which is ranged < e < 1. Here we assume that e is 
a constant for the simplification of our argument, though the actual coefficient of restitution 
depends on the impact velocity [22,23]. The effects of the impact velocity dependence of e 
to macroscopic hydrodynamics can be seen in ref. [24]. It should be noted that the operators 
h and satisfy 

6MG,g) = MG,6g) (8) 



for an arbitrary function h, where 6 is 6 or b~^, and G = (v + Vi)/2. We indicate that our 
model is almost the same as that for inelastic hard-core particles. The difference appears 
through the simplification of the coUisional integral: The term 0(g • 6-)(g • a) for hard- 



core model with the Heviside function Q{x) is replaced by x-^To/m with a constant x- As 
mentioned in introduction this simplification is justified for the potential obeying with 
the relative distance of colliding particles for e = 1. It should be noted that many of inelastic 
Maxwell models contain the factor ^jT/m instead of ^jTo/m in eq.(5) [16,18-20]. However, 



the choice of \jT/m leads to inconsistent results with the elastic Maxwell model in the limit 
of e ^ 1, e.g. the viscosity becomes proportional to \/T. Our choice ensures the consistent 
result with the elastic Maxwell model. Namely, the viscosity is proportional to T. Our 
model may be interpreted as follows. Particles feel the long-range repulsive interaction by 
the potential obeying r~^, and the energy of colliding particles loses by eq.(6) at the instance 
of changing the direction of the relative velocity. 

It is easily verified that the loss of kinetic energy in each coUision is 

l-e^ 



-m(g-<T) . 



(9) 



A useful identity for an arbitrary function for an arbitrary function h is 

j dvh{v)J[f, f]^x'j'J'^J dv, J dvfh J da{g ■a){b- l)h{v), 



(10) 



where / and /i represent /(r, v, t) and /(r, Vi, i), respectively. Prom this identity we directly 
obtain the relations 
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j dVi j dVslVi - V2|V(Vl)/(V2). 



(12) 



The balance equations for hydrodynamic variables are 
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Dtn + nV • u = 0, 



(13) 



DtUi + {mn)-^VjPij = 0. (14) 

and 

AT + ^{PijVjUi + V • q) + TC = 0, (15) 
where Dt = dt + u ■ V , and 

C[/] = (l-e^)^(^[/,/]. (16) 

The pressure tensor P^j and the heat flux q are respectively defined by 

Pij = nT5,j + J dVAi(V)/(r, v, t) (17) 

d^l dVS(V)f(v,v,t), (18) 

where 

D,j{\)^m{V^j-^V%), (19) 

S(V) = (^mV^ - ^r) V. (20) 

Once we know the solution of the inelastic Boltzmann equation, we can obtain the complete 
information of hydrodynamics. It is, however, impossible to get the complete solution of 
the Boltzmann equation which is a nonlinear and a differential integral equation. Thus, we 
need a systematic perturbative scheme to obtain the solution. 

One of the systematic methods to obtain an approximate solution is the Chapman- 
Enskog method [3] . This method is regarded as a standard one in the gas kinetic theory of 
elastically interacting particles. The method assumes a solution of the Boltzmann equation 
to the form 

/(r,v,i) = /[v|n,u,r]. (21) 
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This means that space and time dependences appear through hydrodynamic variables 
n, u, T. We also assume that the solution is nearly homogeneous and exists the small spatial 
variations. Thus, we assume that the distribution function is represented by a series 

/ = f^^ + + eV^^) + • • • , (22) 

where e is a formal expansion parameter and set to be unity after the calculation. The 
expansion parameter is assumed to be balanced with the spatial gradients. In addition, as 
usual in the Chapman- Enskog method for elastically interacting particles, the time derivative 
is also expanded as dt = dt^^^ + edt^^^ + ■ ■ ■. We note dt^^^ = for elastic cases. 

To remove the ambiguity of the distribution function, we impose the solvability conditions 
as usual 

n{r,t)= /(iv/W(r,v,t), (23) 
,t)n{r,t) = I dvv/(°)(r,v,t), (24) 
n(r, t)T(r, t) ^ j dv^mV (r, v, t) . (25) 



3 
2 



Prom these conditions, perturbative distribution functions should be orthogonal to f^^\ 
v/(o) and 



B. The basic solution 

The most difficult part of the Chapman-Enskog method for inelastic particles is to obtain 
the zeroth order solution of the Boltzmann equation 

a^(o)y(o)^ j[^(o)^^(o)]_ (26) 

As mentioned in Introduction, the solution of (26) has been obtained in these days. [18-20] 
In this subsection, we summarize the parts of their results which will be needed for our 
analysis. 

Assuming the scaling form 
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/(v, t) = nvo{t)-'f{c), c = y/vo{t) (27) 
with vo(i) = ^j2T/m, eq.(26) for the inelastic Maxwell model becomes 

--fx^'{x) + ^{x) = ( $(a;e+)$(xe_); cq = -7eo(^), (28) 

J IT 

where = (l/47r) / da and 

e+ = f{k ■ arf, e_ = 1 - ^(fc • (t)^ (29) 

with k = k/A;,p = (l + e)/2 and 2; = l-{l-pf. Here $(eo(t)a;) = 0(a;,i) with eo(t) = 'yo(t)^ 
and <f{k,t) — (f){k'^ /A^t) is the Fourier transform of f^^^n. The solution of (28) including 7 
can be obtained by the combination of the moment expansion and picking up the singularity 
(Appendix A). The result may be written as 



$(a;) = V(-a;)'^^->la;" + ---, (30) 
^0 ^■ 

where a is an non-integer and [a] is the largest integer less than a. The moment //„ is defined 

by 

=< c?^ > I (3/2)„ , (a)„ = r(a + n)/r(a) (31) 

with the Gamma function V{x) and < (?'^ >= J dcc^"/(c). The explicit results for a and 
arc presented in Appendix A. 

The macroscopic equations at the zeroth order are given by 

dt^'^n = 0, dt^^'^u^O, T-%^'^T ^ -C'-''\ (32) 

where the cooling rate ^^^^ is 

C^'^ = {l-e')^u[f^'\f% (33) 

Substituting (32) into (26) we obtain 



The second expression in (34) is based on the assumption that f^^^ is only a function of the 
velocity through the scaled velocity c. 

Here we should stress that f^^^ can be evaluated in the inelastic Maxwell model. The 
solution is determined from eq.(28). It should be noted that the scaling function /(c) defined 
in (A3) has a tail obeying c~^"~^, where the tail is determined in eq.(A7). 

C. The determination of the transport coefficients by the Chapman-Enskog method 

The solution f^^^ is isotropic so that the zeroth order pressure and the heat flux are given 

by 

Pi/'^^pSij, q(°)=0 (35) 

where p — nT is the hydrostatic pressure. 

The first order equation of the Boltzmann equation becomes 

(a,(°) + L)/(^) = -{dP + V ■ V)/(°) = -(A^^) + V • V)/(°) (36) 

with D/^^ = 9/^^ + u • V. Here the hnear operator L in (36) is defined by 

i,j(i)^_j[j(o)^y(i)]_j[y(i)^y(o)]_ (37) 

It is easy to verify that the zero eigenf unctions of (37) are not directly related to the coUisional 
invariants, i.e., /*^°) and v/*^°) are not zero eigenf unctions, but they are zero eigenfunctions 
of This causes significant differences in the perturbation method for systems of inelastic 
particles from those of elastic particles [3,15]. To recover the standard procedure, we need 
to restrict our interest to the case near e = 1. [12] 
Hydrodynamic equations at the first order give 

2T 

DPn= -nV -n, D^n = -{mn)-^Vp, Dt^^^T = V • u. (38) 

3 

Here we have used C^^^ = from the symmetry consideration of variables as in the case of 
the hard-core model. Therefore eq.(36) becomes 
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(dt^''^ + L)/(i) = /(°)(V • u - V • Vlnn) + (9^/^°)) (^^V • u - V • Vt) 

+ (J^/'^) [-(mnr'V,p + Y-Vul (39) 
where Ui is the i-th component of u. The solution of this equation is assumed to be 

/(I) = A(V) • VlnT + B(V) • Vlnn + Cij{V)\/iUj. (40) 

Substituting this into (39), we can determine the functions A,B and Cij. The details of 
calculation is given in Appendix B. The pressure tensor and the heat flux become 

Pij^^^ = -vi'^iUj + VjUi - ^(5i, V • u) (41) 
q(i) = -kVT - i^Vn, (42) 

where r] and k are the shear viscosity and the thermal conductivity, respectively. The other 
transport coefficient /i appears only is granular gases. 

The calculation of the transport coefficients appeared in (41) has been presented in 
Appendices B and C. They are given by 

r;* = ^ = [V-Cr' (43) 
Vo 

, K 2(l + c*(e)) 

« = — = -w? (44) 

«o 3(V-C*) 

Here 770 and Kq are the elastic values of the shear viscosity and the heat conductivity, re- 
spectively. Their values are 



T /2m 

-0 = (47) 

with the numerical coefficient — 1.3700 [25]. Note that the above expression can be 
obtained in terms of the exact perturbative calculation. The collision frequency i^o is defined 
hy p/rjo and its value is 
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uo^3A^na\f^. (48) 
V m 

The dimensionless functions appear in eqs.(43)-(45) are given by 

C(e)^^ = (l-e^)^^^^^^^, (49) 



c*{e) 
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2(/X2-l), (50) 



^ , ^ /dVA,(v)i.Q,(v) .g^. 

" z/o/t^VA,(V)Q,(V)' 
^ , ^ J dVSiV) ■ LAjV) 
" uo J dYSiV) ■ AiVy ^ ^ 

^ /.VS(V).LB(V)_ 3 
i/o/ciVS(V)-S(V) ^ ^ 

To obtain explicit forms of the transport coefficients we need to have the expressions for 
A, B and Cij. The advantage of the inelastic Maxwell model is that we can calculate 
them exactly within this perturbation scheme. From the solvabihty conditions (23)- (25), 
the leading terms of the expression are A oc B oc /^°)S(V) and Cij oc f^°\V)Dij(V). Then 
the final expressions are given by 

C*W = ^(l-e^), (54) 

= JHii^, (55) 

^ ^ 3e2-6e-5' ^ ^ 

. (1±£M^. (56) 

, _ , _ (1 + e)[7A^3(19 - lie) - 20/i2(7 - 3e) + 45 - 5e] 
'^'^ -""^ ~ 24(7/^3 - 10/.2 + 5) ' ^^^^ 

where the explicit expressions of 112 and 113 are presented in (AlO) and (All), respectively. 

Thus, the shear viscosity is given by 

The expressions of k* and /i* are respectively given by 

- ^(^^^ - ^) (59) 
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and 



5(l-e^)(2/X2-l) _ ■ 
12z/ * - 5fl - e2) 



(60) 



where ly,^* and i/^* are given in (57). The behavior of them are given by Figs.l and 2. It is 
easy to verify k* and //* tend to 1 and as e — > 1, respectively. Although they include //a 
which diverges at e = Cc — 0.145123, k* does not have any singularities around the critical 
e, though the value for e < Cc does not have any physical meaning. 

The corresponding first order distribution function is also obtainable. (See Appendix D 
for the derivation). The result is summarized as 



/(i)(V) 



1 



Arrr TlT 



/<»>(V'). (61) 



We are also interested in the stability of uniformly cooling state. For this purpose, we 
need to know the form of As was discussed by Brey et al. [13], the important terms for 
the linear stability analysis is only two terms : 



.(2) 



(62) 



Calculation of and ^2 is possible, as presented in Appendix E. The results are summarized 
as follows: Let us introduce the dimensionless Burnett transport coefficients 



Here Ci* and C2* satisfy 



> * _ 3p _ 3n2 



(63) 



(64) 



where 



with 



, 5,^ 2^/105 15,3 , 9 
= Y2^ ^ V"8~'^' + T^2 + ''^'^^ + 4"" + ^) ' 



5(3/X2 - 7/X3) 

C( = ; P 

2(5//2-3) 



15(/X2^ - 7/X3) 
4(5//2-3) ' 



(65) 



(66) 
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and 



Ct* = — -, , Cn* = — (67) 

Here ct^^^ and Cj^^^ are introduced in (Ell). Prom the calculation presented in Appendix 
E, we obtain ct* and c„* as 

Sk* (63^14 -35^3) _ 16 * 
15(7At3-10A'2+5) 3 f " 

{^e - 3C*) (63/^4 + 14aAt3 + 4/?/X2) - f /«*/X2 

and 



CT* = 777^ o^.n/^o"!!' ! .1! .. ! Aa.\ W-rr (68) 



8m* (63^4-35/^3) _ 16 * 
^ * ^ 15(7At3-10At2+5) 3 

" (i/^*-2C)(63/X4 + 14a/X3 + 4/?//2)-f//*/X2 ^ ^ 

Substituting these results and (65) into (64) we obtain the explicit form of Ci* and ^2*- It 
should be noted that i/^* is given by (E24). 

III. HYDRODYNAMIC EQUATIONS AND THEIR STABILITIES 

The results obtained in the previous section give a closed set of hydrodynamic equations 
for n, u and T at Navier-Stokes order: 

Dtn + nV ■vl = {) (70) 

DtUi + { ViP - V,- [r]{ViUj + VjUi - V • u)] | = (71) 

3n 2 
—DfT + pV ■ u - ViUj[r]{ViUj + VjUi - -SijV ■ u)] 

-V • [kVT + iiVn] = -rC^°^ - TC'-^l (72) 

Here we collect terms up to the second order and set to be e = 1. Since we restrict our 
interest to the linear stability analysis of homogeneous cooling state, as mentioned in the 
previous section, C*^^) can be replaced by Cl*^^^ as in (62). 

For the linear stability analysis, let 5y = y — yn he the deviation from the value in 
a homogeneous state, where y and y^ are a hydrodynamic variable and its homogeneous 
value, respectively. A set of Fourier transformed dimensionless variables are defined by 
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5yk{T) = j d^expHk.^]5y(er), (73) 

^k THT^TT^^, Wk = ./;^^^5Uk T , pkT =— 74 

J-H{r) y iH(r) rii/ 

where the subscript H denotes the quantity at the homogeneous cooUng state. ^ and r are 
the dimensionless space and time variables, 

T^\f dt'Mt') = lA^n^a^l^t, ^ = &, (75) 
2 Jo 2 \ m 2 ^ iH 

Ft' 

with uh = SAjtiUhu'^x — • In terms of these variables the linearized hydrodynamic equations 



are 



drPk + ikwkw = 0, (76) 

{dr -C + lv*k')wk\\ + tk{e^ + Pk) = (77) 

(dr-C + lv*k^)y^k±^0 (78) 

and 

[dr + ^(«* - Ci*)k']e^ + [2C + ^(/^* - C2*)k'W + likWkw = 0. (79) 

Here Wk\\ and w^^ denote the longitudinal and the transversal component of Wk defined in 
eq.(74), respectively. 

Equation (78) is decoupled from the rest and the solution is given by 

Wfe±(T) = Wjk_L(0) exp[s±T] (80) 

where 

s±^C- \ri*k\ (81) 

This identified the degenerated shear modes. The remaining eigenmodes have the form 
exp[s„T] for n = 1, 2, 3, where s„ are the solutions of 
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(82) 



The dispersion relation is summarized in Fig. 3. This result indicates that two of real parts 
of Sn are degenerated in all regions. In addition, the shear mode does not have the largest 
eigenvalue in the unstable region. 

The hnear hydrodynamics discussed here contains the characteristic wavelength; 



initial long wavelength perturbations of homogeneous cooling states grow exponentially. To 
discuss the spatial structure after the instability of the homogeneous state, we need a proper 
set of nonlinear equations. Although we do not know what the proper nonlinear equation is, 
we may expect that the time-dependent Ginzburg-Landau type equation may be a candidate 
to characterize a nonlinear region [26]. 

In any case, this result suggests that the system of granular gases does not have any 
entropy. So the expectation that the Tsallis entropy [27] can be used in granular systems is 
hopeless. 



In this paper, we have derived hydrodynamic equations based on the Chapman- Enskog 
method. The analysis presented here is so systematic and straightforward that we can discuss 
problems in hydrodynamics of inelastic particles in more general situations, e.g. how we can 
apply hydrodynamics to granular systems. 

First of all, our result indicates that a transport coefficient /i diverges at ec as a result of 
divergence of /i^. This divergences is closely related to the high energy tail of VDF obeying 




(83) 



in which the uniform spatial structure larger than (A;±'^) 



^ and {k\\'^) ^ is unstable. Therefore, 



IV. DISCUSSION 
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a power law in the homogeneous coohng state. Our analysis contains the information of the 
distribution function completely which cannot be achieved for systems of inelastic hard-core 
particles. 

Second, we can consider a driven system by adding the Fokker-Planck operator [28,29] 

df 

^ + V • V/ = J[f, f] + Lppf, (84) 

where 

Lfp^^o^-[-v+—-^]. (85) 
av m ov 

This model is physical, because the particles feel fluid drag —70V and the thermal activation 
by the heat bath Tg. This system reaches a steady state, and as the scaling solution in 
homogeneous cooling states is equivalent to a steady solution in the Gaussian thermostat 
system in eqs.(84) and (85) at Tb = [30]. It is obvious that the system at finite Tb still 
has similar properties to those in cooling systems. Although it has been recognized that 
patterns of free cooling systems [31] are similar to those in suspension of fluidized beds 
[32], no connection has been discussed systematically. In fluidized beds, there is the long 
ranged hydrodynamic interactions among particles [21], which means hard-core model in 
the coUisional integral may not be appropriate for the system with flow. In any case, the 
analysis presented here is not only limited to mathematical interest which can be solved 
in terms of the exact perturbation but also may cause physical interest in applications to 
granular particles in fluid flows. We note that the high energy tails obeying a power law in 
VDF have been reported in some other papers of granular systems [33,34], but the model 
in one of them [33] may not be appropriate as recognized by the authors themselves of the 
paper, and another results [34] could not be reproduced by any other groups. 

Third, one can ask if the granular temperature introduced in (3) is an actual hydrody- 
namic variable. For e — 1, mV'^f^^^ oc c^f{c) is the zero eigenfunction of L. As a result, T 
for a spatial homogeneous state does not have any relaxation mechanism. In our dissipative 
system, T is not a true hydrodynamic variable but a quasi-hydrodynamic variable in strict 
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sense. In other words we still do not have clear picture why we can assume the solvability 
condition for T as in the case of e = 1. We, at least, need to show that the separation 
of the eigenvalues of L of c^f{c) from other modes. At present, though we do not have 
any proof of this requirement, we may have rough picture to support this. Let us consider 
analytic functions which can be represented by polynomials of c. What we need to show is 
the eigenvalue of i?„(c)/(c) with n th. order polynomial is much larger than that of c^ f{c). 
This seems to be true, because oc / dcRn{c)L[Rn{c) f{c)\ ~ We remember that there 
is an inequality > [in-i and its equality is achieved at e = 1. Since higher order mo- 
ments diverge in our model. Even when /x„ is finite, fin with large n is much larger than 
A*2 > A*i = 1- Prom this simple consideration, we can expect the separation of eigenvalues 
between energy and the others. As a result we may justify to assume that the temperature 
is a hydrodynamic variable. 

If we believe that non-Gaussian properties or the violations of the detailed balance are 
essential to steady states or scaling region of cooling states in general dissipative systems 
besides granular gases, we need investigations which is not based on expansions around 
Gaussian distribution function. In this sense, our systematic analysis presented here is 
a good example to demonstrate how non-Gaussian nature affects macroscopic behaviors 
of the system. We hope that the analysis of this 'solvable' model gives some insights to 
understand macroscopic behaviors in general situations of non-Gaussian systems or locally 
nonequilibrium systems. 

V. CONCLUDING REMARKS 

We have derived hydrodynamic equations for the inelastic Maxwell model based on a 
systematic Chapman- Enskog method. We have determined all of the transport coefficients 
77, n and n appear in Navier-Stokes order as a function of the restitution coefficient e. They 
cannot be defined for small e as a result of the high energy tail of the velocity distribution 
function obeying a power law. We also determine the dispersion relations for the linearized 
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equation around the homogeneous coohng states. Finally, through the analysis in this pa- 
per, we clarify the limitation of conventional gas kinetics method, i.e. the hydrodynamics 
corresponding to the Navier-Stokes equation exists only in e > Cc- 

Note added after submission: Immediately after my submission, I have recieved preprint 
whose subject is closely related to that of this paper [35] . Santos [35] adopts the inelastic 
Maxwell model with a prefactor ^T/m instead of ^To/m. He also discusses steady states 
by adding thermostats. Although one can see some equivalent results between mine and 
his, differences exist in (i) the model is different, and (ii) linear stability of a homogeneous 
cooling state is discussed in this paper. 

The author would like to thank N. Mitarai for fruitful discussion. He also appreciate 
useful comments by S. Sasa and H-D. Kim. He expresses his sincere gratitude to Prof. 
A. Santos who let him know ref. [35]. This work is partially supported by the Hosokawa 
Powder Technology Foundation, and the Inamori Foundation. This paper is dedicated to 
the memory of Daniel C. Hong who has passed away in July, 2002. The author got an idea 
of this research during his stay at Lehigh University through discussion with D. C. Hong. 
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FIG. 2. /X* as a function of e. 
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FIG. 3. The dispersion relation as a function of dimensionless k at e = 0.9. s{k) denotes the 
dispersion of the shear mode. 

APPENDIX A: THE DETERMINATION OF THE BASIC SOLUTION 

In this Appendix we summarize the method of determination of f^^^ obtained by some 
authors in these days. Now we introduce 



where /* and f*^ arc, respectively, the precollisional / and /i, i.e. /* = /(v — ^^(g ■ o')a) 
and /*! — f{vi + ^(g • Let us assume the scahng form of VDF as 
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(Al) 



Thus, eq.(4) becomes 




(A2) 



/ = ^o(r)-V>M(r)) (A3) 
with Vo{t) — ^j2T/m. The scahng function /(c) satisfies the normahzation 

Jdcf{c)^l, Jdcc'fic)^^. (A4) 



With the aid of (A4) substituting (A3) into (A2) we obtain 

vo^-ivo, 7^ -(c/^)) = /(/,/), (A5) 

where 7 is the separation constant. RecaUing (Al) the time evolution of vo{t) — vo{t) is 
obtained as 



Vo{t) =Vo{0)e-^^ = Vo{0)exp 



2 To. 
-Xna \ —t 
V m 



(A6) 



Since the coUisional integral in the inelastic Maxwell model is independent of the relative 
speed, Vq does not have an algebraic decay but has the exponential decay (A6). 

Bobylev et al. [16] indicate that the Fourier transform of this model becomes an easy 
equation to discuss the behavior. Introducing the Fourier transform (p{k,t) = 0(|/c^,t) of 
f{v,t) and its scaling form 4>{x,t) = ^{eo{t)x) with eo{t) = t'o(t)^ satisfies (28). 

This model has the singularity near x — oi ^{x) in eq. (30) and the moments with 
n > a diverge. To determine a we use 

l = h{l-p)= [[l-e+-e_], a7= /[l-e+"-e_«] (A7) 

o J a J & 

with /- = (l/47r) / do'. From eq.(A7) the exponent can be determined. This singular term 
reflects on the tail of VDF obeying a power law /(c) ~ c~^"~^. 

The moment //„ in (31) which is introduced in eq.(31) satisfies the iterative equation 

-j^ n—l 

lj,„ = — J2 H{k, n - k)iin-klJ'k, (A8) 

7n k=l 

where /ii — 1 and ; 
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H{m, n) — 



( \ 

m + n 



m 



J 



(32mP''^F{-n, m + 1/2, m + 3/2; z) 



m + n 


n 

E 




[ m ) 


k=0 





(A9) 



with the Gaussian Hypergeometric function F{a,b,c;z) = J2n ^ and ~ 

(1/2)^/(3/2)^ = 1/(1 + 2m). From this expression, we can write the exphcit forms of 
Ii2, A*3 and 114 as 

11 - 6e + 3e2 



5 + 6e — 3e^ ' 

115 - 54e + 66e2 - 30e=^ + 15e^ 
33 - 242e + 1026^ - lOe^ + Se^ ' 



(AlO) 



(All) 



and 



where 



/.4 = -[A^3(i^(l, 3) + H{3, 1)) + fX2'H{2, 2)] 

74 



(A12) 



and 



(A13) 



^(3,l) = 4/[^-^], 



(AM) 



1 2 ^2 



/7(2,2) = 6p^(---z+-) 



(A15) 



4 6 ^ 4z^ 8 ^ 

'14, = -z z H pil — p). 

'35 7 9 9 3^^-^^ 



(A16) 



We note that fin are monotonically decreasing functions of e, and are reduced to = 1 
at e = 1. 112 exists for all e but 112 and /X3 diverge at 0.145123.. and 0.38386.., respectively. 
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APPENDIX B: CHAPMAN-ENSKOG SOLUTION 



The Chapman-Enskog method to the Boltzmann equation is an estabhshed method to 
obtain an asymptotically correct solution as a series of spatial gradients [3,15]. In this 
Appendix we explain the method to obtain the Chapman-Enskog solution. It should be 
noted that this Appendix is parallel to Appendix A of Brey et al. [13] for hard-core inelastic 
gases. The differences appear in the evaluation of the viscosity and the thermal conductivity 
for elastic Maxwell model. 

The first order equation of f^^^ may be written as 

(9/°) + L)/(^) = A- Vlnr + B- Vlnn + Qj^jUi. (Bl) 
Substituting (34) and (35) into (39), and comparison it with (Bl) leads to 

A(VKT) = |A.(v/(o.)-l^/(o.. (B2) 

B(V|«, T) = -V/l»> - -^f'"' (B3) 

and 

QAy\n,T) = -§^iy,f'^) - ■ (V/W). (B4) 

The solution of eq.(Bl) is assumed to be the form of (40). We note again that C^^^ = 
directly from the symmetric considerations, since C^^^ is a scalar. We also indicate the 
relation ^t^^^VlnT = -VC^°^ = -C^^^Vlnn. 

From (40) and (Bl) with the help of the last equation and (32) we reach 

(-C^^^T^T + l)a = A (B5) 

{-C^^^Tdr + l) B = B + C^^^A (B6) 

(_^io)Tdr + L)Cij = Ci^. (B7) 
24 



Let us calculate the viscosity at first. The pressure tensor at the first order is written as 



3 

After a long calculation we can show the relation 



2 

-vi^jUi + ViUj - -6ijW ■ u). (B8) 



10 

Introducing 

_ /rfVA,(V)LC',,(V 



V--^ I dVD,,{Y)Q,{V). (B9) 



" JdYD,j{\)C,j{y) 

eq.(B7) is reduced to 



(BIO) 



(-C(^)TdT + i^,)v --^J dyD,j(V)aj{Y). (Bll) 
This equation can be solved as 

^^-ToK^/^^^^^^ 

^ ^ dVmV^f'^ = (B12) 



In the elastic Maxwell model in the dilute gas, the viscosity is given by [1,3,25] 



1 2m T 1 2m p 

^°=3V^a;;^-3V^^^ ^^''^ 

and there is a relation between r]Q and the thermal conductivity Kq as 

where — 1.3700 is a constant, and To is the strength of the repulsive potential. Let us 
introduce a characteristic collision frequency defined by 

I/O = — = SA^na^J—. (B15) 
r)o \ rn 

Thus, (B12) becomes 
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V _ 1 
Vo - C* 

where Vj^* — V'n/^o and C* = C^^V^o- 
The heat flux of this order is 

q(i) = -kVT - nVn 

with the transport coefiicients 

K = -^/ <ivs(v) . A(V), 

Similar to the viscosity, k and obey 

.^^^j^/dVS(V).A(V). 



1 

1^ 



f„ - 2C<°> 
where 

/ f/VS(V) • I A(V) / r/VS(V) • IB(V) 



J (iVS(V) • A(V) J (iVS(V) • B{Y) 

Use of the formula for A and B the further simpliflcations are given 

-/rfVS(V).A(V)=-— [l + c-(e)l 
i-/,VS(V).S(V)^-5Z|M 

where c*(e) is given by 

c*(e) = 2(//2 - 1) 

with the moment /i2 introduced in (AlO). 
Prom these results k and // become 

^ ^ ^ 2 (l + c*(e)) 
3 (i/k* - C*) ' 

* = = C*^/^o + c*(e)/3 
where i/^* = i^k/^^o and i/^* = i^^^/i^o- 
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APPENDIX C: EVALUATION OF TRANSPORT COEFFICIENTS 



To evaluate dimensionless transport coefficients (*, ^ and we need to know 
the details form of distribution function discussed in Appendix A. We should note that we 
cannot separate integrals of V and Vi by G and g, since basic VDF for homogeneous states 
is not the Maxwellian. 

First, we evaluate C*- Prom (49) has the form 



,,2u;|/<«l,/<»)] 



Prom (12) and (A3) with (Al) we obtain 

a;(/,/) = ^^'"^^^ / dci / cic^lc- Ci|V>)/(ci). (C2) 
Here we use / dcc/(c) = 0. Substituting this into (CI) with the help of (A4) we obtain 

6Vq V ra 

The other transport coefficients ^ v^* and i/^* are evaluated as follows. To lowest-order 
velocity dependence is 

A(V)cx/WS(V), S(V)«/(°)S(V), A,(V)oc/(°)A,(V). (C4) 

Prom (B21) and the definition of v* with i — r], k, jj, we obtain 

JdYD,,{Y)L[f('^D,,iY)] 
uoJd\p)Dij{\)Dij{\) 

^ /c^VA.(V)L[/WA,(V)] (C5) 



10nT2(l + c*(e)/2) 



and 



Here we have used 



, _ /rfVS(V)-L[/WS(V)] 
^/o/rfV/(o)S(V).S(V) 
4m/(iVS(V) ■L[/WS(V)] 
ISnTVo (7/^3 - 10/^2 + 5) ' 
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(C6) 



dV/(°)S(V) ■ S(V) = ^— (7/X3 - 10/X2 + 5) . (C7) 

Prom the linear collision operator defined in eq.(37), the integral of the form / YL[f^'^^X] 
can be transformed as 



2 1^0 

V m 

If. 



dvr(v)L[X(v)/(°)] 

j dw,j dw J daY{w){e-'h-' - l)/(°)(v)/(°)(vi)(X(v) +X(vi)) 

= / rfvi |rfv| da/W(v)/(°)(vOX(vi)(6- l)(F(v) + F(vO) 

= I civiX(vi)Lt[r(v)/(°)]. (C8) 

The evaluation of the integrals of eqs.(C5) and (C6) is now possible from the straight- 
forward calculation. Let us calculate v^* at first. Noting Dij(V)Sij — 0, eq.(C5) can be 
replaced by 

d^r J dvi/('^)(V)/('^)(Vi)A,(Vi) I d&{b-l)iViVj + VuV,j) 

(C9) 



~ 10nz/oT2(l + c*(e)/2) 



Following the textbook by Chapman and Cowling [3] , the solid angle integral over a can be 
performed. Noting (6 - l){ViVj + VuVij) = ^gkC^k{gm^i(^j{'^ + e) - Qi^j - Qj^i} and 

An 



/47r 
d&{b - l){ViVj + VuV.j) = — (1 + e)(e - 4) A,(Vi)y,^,- (CIO) 



we obtain 



27ra2xV^(l + e)(4-e) 



JdvJ dVi/(°)(V)/(°)(Vi)A,(Vi)^,^,-. (Cll) 



" TSnTVo (l + c*(e)/2) 
The result of the integration over V and Vi leads to 

In the limit of e = 1, i/^* = 1 should be recovered. Thus, parameters x should satisfy 

X = = 0.8176582 ••• , (CIS) 
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which is deviated from 1 a httle. Thus, with the aid of (48) the final expression becomes eq. 
(56). This leads to the expression of rj* as in eq.(58) 

On the other hand, the evaluation of u^* and i/^* are as follows. From the definition of 
S(V) in eq.(19) and the coUisional operator b in (7) / dVS(V)L[/(°)S(V)] contains the solid 
angle integral 

J rf(T(6-l)[S(V) + S(Vi)] 

= - + ^) / ^^[(1 - - (1 + e)(g • ^)^(G . a)a- 

+(g-a){(G-(T)g + (G-g)6-}] 
= ^(1 + e)[(l + e)g(/G + 2(g • G)g) - |(1 - e)/G - |g(G • g)] 
- '''"^^ + '^{(3-7e)52G + 4(4-e)(g-G)g}, (C14) 
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where G = (V + Vi)/2. Thus, the numerator of (C6) becomes 

J dVS(V) • L[/WS(V)] = ^^^j + '^ a^^ /rfv/ rfVi/(°)(V)/(°)(VO 

S(Vi) • {(3 - 7e)/G + 4(4 - e)(g • G)g}. (C15) 

Here we can rewrite 

/ dVS(V)L[/WS(V)] = ^I^^^Mt^I + e)[^^h + 2(4 - e)h] (C16) 
where /i and /2 are respectively given by 

7i = ycic y'dci/(c)/(ci)(ci^-5/2)(c-Ci + ci^)|c-Ci|^ 

105 15 15 
- ^/^3 - + y (C17) 

and 

h^-jdcj dcj{c)f{c,){c,^ - 5/2) (c^ - cl)c^ 

105 45 , 

= — /X3- 15/^2 + y. (CIS) 

Substituting (C17) and (CIS) into (C16) we obtain 
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/ (iVS(V)L[/(°)S(V)] = + e)[7Ml9 - lie) - 20/^2(7 - 3e) + 45 - 5e]. 

7 4m V JT^' 

(C19) 

Substituting this into (C6) we reach eq.(57) . 

APPENDIX D: THE DISTRIBUTION FUNCTION AT THE FIRST ORDER 

The exphcit form of the distribution function at the first order is needed for the discussion 
of hnear stabihty of homogeneous coohng states. Let us assume the form 

/(^) = [ctS(V) • VlnT + c„S(V) • Vlnn + c„A,(V)Vi«,]/(°)(V). (Dl) 

Here recaUing eqs.(40),(50) and (B9) we can rewrite 

"^^'Tol ^VA,(V)Q,(V) 

X5 J 

= -nT^c^H2- (D2) 

From (B18), similarly, we obtain 

= -^/rfVS(V)-S(V)/W 

^''^''^'[7/.3-10/X2 + 5] (D3) 



4m 



and 



= ~ I (iVS(V) . S(V)/<°1 



^^-^[7/i3-10/i2 + 5] (D4) 



4m 



Therefore we obtain /(^) as 

/a)(v): 



S(V) • {kVT + iiVn) + ^Dij{y)ViUj 



5(7//3 - 10//2 + 5) //2 



f<^\V). (D5) 



30 



APPENDIX E: DETERMINATION OF C^^^ 

In this section, let us determine C^^^ which is needed to discuss the hydrodynamic stabihty 
of homogeneous states. The second order correction C*^^^ now becomes 

C^^^ = (1 - e^)^ [^[f^'\ f^'^] + Mf'\ f'^]] (El) 
Now let (l^"^^ be the hnear part of which satisfies 

CL^'^ = i^-e')^u;[f'\f% (E2) 

The second order equation now becomes 

idt'-'^ + L)f^'^ = - (A^^) + V • V)/« + J[/W, /«]. (E3) 

Here the contribution from J[f^^\ /*^^^] can be neglected, because this term creates nonlinear 
terms of hydrodynamic variables. We also note 

dt^'^^n = 
mndt^'^^u = ^V(V • u) + rjV'^u 

^^(2)^ ^ _T({2) ^ A(«;V2r + iiV^n) (E4) 
3n 

for linearized equations. Taking into account the above argument [3] and (38) we obtain 

= -|^(«v2r + /xv^n)aT/(°) 

4mS(V)V : (kVVT + iiVVn) 
^ 5nT%7n3 - 10//2 + 5) 

+^llv(V.u) + V^u]. A;(o) 
nm 3 o\ 
Am 

~5nr3(7/X3- 10/X2 + 5) V3^ 
-^^^A,(V)[(nm)-V,V,p - V : V(V,m,)]/(°). (E5) 

Further simplification is possible by noting that (^^^ is a scalar so that any contributions to 

/l*-^-* that are vector or traceless functions does not contribute. Let Sfi^^^ be the residual 
part oifL^^l Thus, we can simplify 
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'-Tk + ni?j S(V) • (VV • u) 



. 4m /W(V)S(V)-V 



Let us assume 



15nT^ {TfjLs - IO/X2 + 5) SnT 



5fiP = M(r, V) V^T + N{T, \)V'n 



TdTf^^\K\7^T + jiV^n). (E6) 



(E7) 



and taking into account (62) we obtain 



(0) 



'n. 



n 



(E8) 



Substituting this into (E7) eq.(E6) becomes 



Am /(o)S(V)-V 2 



and 



15nr3 7/^3 - IO/X2 + 5 3nr 
4m /(o)S(V)-V 2 



(E9) 



(ElO) 



The scalar functions M and A^" are orthogonal to 1, v and from the solvabihty conditions 
(23)-(25) , i.e. not to contribute the transport coefficients determined in Appendix C. Thus, 
we can assume 



(Ell) 



with the form 



Q(c) = + ac^ + (3. 



(E12) 



Prom the orthogonal conditions, a and (3 can be determined as (66). Thus, we obtain (65). 

Equations to derive c^^^^ and Cn^"^^ are obtained from eqs.(E9) and (ElO). At first they 
are rewritten as 



8/c2(c2-5/2) d 



5 (7/X3 - 10/^2 + 5) dc 



(E13) 
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and 



n 



3nTi/o 



8/V(c^ - 5/2)) ^ 
5(7//3 - 10//2 + 5) dc 



(c/) 



Multiplying them by and integrate them all over c we obtain 
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63//4 — 35/^3 



15pi/o(7//3 - 10//2 + 5) 63//4 + 14q;//3 + 4/3//2 



(E14) 



(E15) 



i^o V 2 3py 63/^4 + 14q;//3 + 4(3/12 

jj, 63/X4 — 35/i3 



15pi/o(7/X3 - 10/^2 + 5) 63/^4 + 14q;/X3 + 4/3/^2' 



where 



/ricc^L[Q(c)/(c)] 
i/o/rfcc4g(c)/(c) 

4iy 



(E16) 



(E17) 



Here is defined as Vl^ = / o?cc^L[(5(c)/(c)]/i/o- Introducing dimensionless ct* and c„* by 



Ct = 



Ct 



(2) 



(2) 



eqs.(E15) and (E16) can be rewritten as 



(i/^* - 3C*)(63/i4 + 14q;//3 + 4/3//2) 

X 



8(63/^4 - 35//3)/^* 16/^2 ^ *x 
15(7/^3 - 10/^2 + 5) 3 ' 



and 



(E18) 



(E19) 



{u^* - 2C*)(63At4 + Uans + 4/3/^2) 

8(63/7,1 — 35fi3)fi* ^ 16/i2 



X 



15(7At3 - 10/X2 + 5) 



+ 



(C2* - K*) 



(E20) 
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Eliminating Ci* and C2* from (E19) and (E20) we obtain (64). 

Now, we can evaluate 1^ as in the previous section. The result may be 



W 



8 



(E21) 



where 



7 -^^ 



,15 3 63 7 



(E22) 



and 



3 /105 15 3^ 
^ 2 l^^^ + T"^^ + 2^ 



(E23) 



Thus i/^* is determined by 



W 



15 f /X4 + l^yUa + P1JL2' 



(E24) 
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